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Abstract. A framework is a graph and a map from its vertices to E d (for some d). A 
framework is universally rigid if any framework in any dimension with the same graph 
and edge lengths is a Euclidean image of it. We show that a generic universally rigid 
framework has a positive semi-definite stress matrix of maximal rank. Connelly showed 
that the existence of such a positive semi-definite stress matrix is sufficient for universal 
rigidity, so this provides a characterization of universal rigidity for generic frameworks. We 
also extend our argument to give a new result on the genericity of strict complementarity 
in semidefinite programming. 



1. Introduction 

In this paper we characterize generic frameworks which are universally rigid in d-di- 
mensional Euclidean space. A framework is universally rigid if, modulo Euclidean trans- 
forms, there is no other framework of the same graph in any dimension that has the 
same edge lengths. A series of papers by Connelly [Con82, Con99, ConOl] and later Al- 
fakih [Alf07a, Alf07b] described a sufficient condition for a generic framework to be univer- 
sally rigid. In this paper we show that this is also a necessary condition. 

Universally rigid frameworks are especially relevant when applying semidefinite program- 
ming techniques to graph embedding problems. Suppose some vertices are embedded in E d 
and we are told the distances between some of the pairs of vertices. The graph embedding 
problem is to compute the embedding (up to an unknown Euclidean transformation) from 
the data. This problem is computationally difficult as the graph embeddablity question is in 
general NP-HARD [Sax79], but because of its utility, many heuristics have been attempted. 
One approach is to use semidefinite programming to find an embedding [LLR95], but such 
methods are not able to find the solutions that are specifically d-dimensional, rather than 
embedded in some larger dimensional space. However, when the underlying framework is, in 
fact, universally rigid, then the distance data itself automatically constrains the dimension, 
and therefore the correct answer will be (approximately) found by the semidefinite program. 
The connection between rigidity and semi-definite programming was first explored by So and 
Ye [SY07]. 

We also discuss the more general topic of strict complementarity in semidefinite pro- 
gramming. Strict complementarity is a strong form of duality that is needed for the fast 
convergence of many interior point optimization algorithms. Using our arguments from uni- 
versal rigidity, we show that if the semidefinite program has a sufficiently generic primal 
solution, then it must satisfy strict complementarity. This is in distinction with previous 
results on strict complementarity [AH097, PT01] that require the actual parameters of the 
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program to be generic. In particular, our result applies to programs where the solution is of 
a lower rank than would be found generically. 

1.1. Rigidity definitions. 

Definition 1.1. A graph T is a set of v vertices V(T) and e edges S(T), where S(T) is a set 
of two-element subsets of V(T). We will typically drop the graph T from this notation. A 
configuration p is a mapping from V to E". Let C(V) be the space of configurations. For 
p G C(V) and u G V, let p{u) denote the image of u under p. Let C d (V) denote the space 
of configurations that lie entirely in the E d , contained as the first d dimensions of E 1 '. A 
framework (p, T) is the pair of a graph and a configuration of its vertices. For a given graph T 
the length-squared function £r : C(V) — > R e is the function assigning to each edge of V its 
squared edge length in the framework. That is, the component of £r(p) in the direction of 
an edge {u, w} is \p{u) — p(w) | 2 . 

Definition 1.2. A configuration in C d (V) is proper if it does not lie in any affine subspace 
of K d of dimension less than d. It is generic if its first d coordinates (i.e., the coordinates 
not constrained to be 0) do not satisfy any algebraic equation with rational coefficients. 

Remark 1.3. A generic configuration in C (V) with at least d + 1 vertices is proper. 

Definition 1.4. The configurations p, q in C(V) are congruent if they are related by an 
element of the group of Eucl(w) of rigid motions of E 1 '. 

A framework (p, T) with p G C d (V) is universally rigid if any other configuration in C(V) 
with the same edge lengths under £ r is a configuration congruent to p. 

A graph T is generically universally rigid in E d if any generic framework (p, T) with 
p G C d (V) is universally rigid. 

Essentially, universal rigidity means that the lengths of the edges of p are consistent with 
essentially only one embedding of V in any dimension, up to v. (In general, a configuration 
in any higher dimension can be related by a rigid motion to a configuration in K v .) This is 
stronger than global rigidity, where the lengths fully determine the embedding in the smaller 
space E d . And global rigidity is, in turn, stronger than (local) rigidity, which only rules out 
continuous flexes in E d that preserve edge lengths. 

Remark 1.5. Universal rigidity or closely related notions have also been called dimensional 
rigidity [Alf07a], uniquely localizable [SY07], and super stability [Con99]. 

Remark 1.6. In Definition 1.4, it would be equivalent to require that (p, V) be (locally) rigid 
in K v since by a result of Bezdek and Connelly [BC04] any two frameworks with the same 
edge lengths can be connected by a smooth path in a sufficiently large dimension. 

Remark 1.7. Universal rigidity, unlike local and global rigidity, is not a generic property of 
a graph: for many graphs, some generic frameworks are universally rigid and others are not. 
For instance, an embedding of a 4-cycle in the line E 1 is universally rigid iff one side is long 
in the sense that its length is equal to the sum of the lengths of the others. On the other 
hand, some graphs, such as a simplex, and any trilateration graph [EGW + 04] are generically 
universally rigid in E d . (A <i-trilateration graph has the property that one can order its 
vertices such that the following holds: the first d + 2 vertices are part of a simplex in T 
and each subsequent vertex is adjacent in V to d + 1 previous vertices.) We do not know a 
characterization of generically universally rigid graphs. 

See Figure 1 for examples of embedded graphs with increasing types of rigidity. 
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FIGURE 1. Planar graph embeddings with increasing types of rigidity. From 
left to right, generically locally flexible graph, generically locally rigid graph, 
generically globally rigid graph, generic universally rigid framework, generi- 
cally universally rigid graph. (More precisely, the fourth framework is not 
itself generic, but it and any generic framework close to it are universally 
rigid.) This paper focuses on frameworks of the type of the fourth one. The 
fifth graph is generically universally rigid as it is a 2-trilateration graph. 

1.2. Equilibrium stresses. 

Definition 1.8. An equilibrium stress matrix of a framework (p, T) is a matrix Q indexed 
by V x V so that 

(1) for all u,w G V, Q(u,w) = Q(w,u); 

(2) for all u, w G V with u ^ w and {u, w} G" £, Q(u, w) = 0; 

(3) for all u G V, YlweV w ) = 0> an d 

(4) for all u G V, E«, e v w)p(w) = 0. 

(The last condition is the equilibrium condition.) Let Sr(p) be the linear space of equilibrium 
stress matrices for (p, T). 

Conditions (3) and (4) give us: 

(1.9) VmgV tt{u,w)(p{w)-p(u)) = 0. 

{weV\{u,w}££} 

The kernel of an equilibrium stress matrix Q of a framework (p, T) always contains the 
subspace of spanned by the coordinates of p along each axis and the vector 1 of all 
l's. This corresponds to the fact that any affine image of p satisfies all of the equilibrium 
stresses in Sr{p)- If p is a proper d-dimensional configuration, these kernel vectors span a 
(d + l)-dimensional space, so for such frameworks rankfi < v — d — 1. 

Definition 1.10. We say that the edge directions of (p, T) with p G C d (V) are on a conic 
at infinity if there exists a symmetric d-by-d matrix Q such that for all edges (w, v) of T, we 
have 

\p(u) - p{v)fQ\p{u) -p{v)} = 0, 

where the square brackets mean the projection of a vector in E 1 ' to E d (i.e., dropping the 0's 
at the end of p{u)). 

Remark 1.11. If the edges of (p, T) are not on a conic at infinity in C d (V), then in particular 
p is proper. 

A framework has edges on a conic at infinity iff there is a continuous family of <i-dimensional 
non-orthonormal affine transforms that preserve all of the edge lengths. However, if the 
configuration is not proper, this family of affine transforms might all be congruent to the 
original one. 
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1.3. Results. Connelly in a series of papers has studied the relationship between various 
forms of rigidity and equilibrium stress matrices. In particular he proved the following 
theorem that gives a sufficient condition for a framework to be universally rigid. 

Theorem 1.12 (Connelly). Suppose that T is a graph with d + 2 or more vertices and p is a 
(generic or not) configuration in C d (V). Suppose that there is an equilibrium stress matrix 
Q £ Sr{p) that is positive semidefinite (PSD) and that rankf2 = v — d—1. Also suppose that 
the edge directions of p do not lie on a conic at infinity. Then (p,T) is universally rigid. 

Connelly also proved a lemma that allows him to ignore the conic at infinity issue when p 
is generic. 

Lemma 1.13 (Connelly). Suppose that T is a graph with d + 2 or more vertices and p is a 
generic configuration in C d (V). Suppose that there is an equilibrium stress matrix Q £ Sr(p) 
with rank v — d — 1. Then the edge directions of p do not lie on a conic at infinity. 

Putting these together, we can summarize this as 

Corollary 1.14. Suppose that V is a graph with d + 2 or more vertices and p is a generic 
configuration in C (V). Suppose that there is a PSD equilibrium stress matrix Q £ Sr(p) 
with rankfi = v — d — 1. Then (p, T) is universally rigid. 

The basic ideas for the proof of Theorem 1.12 appear in [Con82] where they were applied 
to show the universal rigidity of Cauchy polygons. It is also described in [Con99] and stated 
precisely in [ConOl, Theorem 2.6]. Lemma 1.13 can be derived from the proof of Theorem 1.3 
in [Con05]. Corollary 1.14 in two dimensions is summarized in Jordan and Szabadka [JS09]. 

A different set of sufficient conditions for the related concept of dimensional rigidity were 
described by Alfakih [Alf07a]. In [Alf07b] he showed that these conditions were equivalent 
to the existence of a PSD equilibrium stress matrix of maximal rank. In [Alf07b] he also 
showed that a generic dimensionally rigid framework (with v > d + 2) must be universally 
rigid; this has the same effect as Lemma 1.13, and thus results in Corollary 1.14. In these 
papers, he also conjectured that for a generic universally rigid framework, a maximal rank 
PSD equilibrium stress matrix must exist. 

In the related context of frameworks with pinned anchor vertices, So and Ye [SY07] showed 
the appropriate analogue of Theorem 1.12 follows from complementarity in semidefinite 
programming (see section 5 below for more on this). 

In this paper our main result is the converse to Corollary 1.14: 

Theorem 1. A universally rigid framework (p, T), with p generic in C d (V) and having d + 2 
or more vertices, has a PSD equilibrium stress matrix with rank v — d—1. 

Remark 1.15. Alfakih has given an example [Alf07a, Example 3.1] showing that Theorem 1 
is false if we drop the assumption that p is generic. For any universally rigid framework 
(generic or not), it's not hard to see that there is a non-zero PSD equilibrium stress matrix. 
(See [Alf09, Theorem 5.1].) The difficulty in Theorem 1 is finding a stress matrix of high 
rank. 

Theorems 1.12 and 1 compare nicely with the situation for global rigidity, where a frame- 
work is generically globally rigid iff it has an equilibrium stress matrix of rank v — d—1 (with 
no PSD constraint). Sufficiency was proved by Connelly [Con05], and necessity was proved 
by the authors and Alex Healy [GHT10]. 
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Question 1.16. For a graph T that is generic globally rigid E d , is there always a generic 
framework in C d (V) that is universally rigid? 

In this paper, we also prove more general versions of Theorem 1 in the general context 
of convex optimization and strict complementarity in semidefmite programming. See Theo- 
rem 2 in Section 4 and Theorem 3 in Section 5. 

2. Algebraic geometry preliminaries 

We start with some preliminaries about semi-algebraic sets from real algebraic geometry, 
somewhat specialized to our particular case. For a general reference, see, for instance, the 
book by Bochnak, Coste, and Roy [BCR98]. 

Definition 2.1. An affine real algebraic set or variety V contained in R n is a subset of M. n 
that is defined by a set of algebraic equations. It is defined over Q if the equations can be 
taken to have rational coefficients. An algebraic set has a dimension dim(V), which we will 
define as the largest t for which there is an open subset of V (in the Euclidean topology) 
homeomorphic to R*. 

Definition 2.2. A semi- algebraic set S is a subset of R n defined by algebraic equalities 
and inequalities; alternatively (by a non-trivial theorem), it is the image of an algebraic set 
(defined only by equalities) under an algebraic map. It is defined over Q if the equalities 
and inequalities have rational coefficients. Like algebraic sets, a semi-algebraic set has a well 
defined (maximal) dimension t. 

We next define genericity in larger generality and give some basic properties. 

Definition 2.3. A point in a (semi-) algebraic set V defined over Q is generic if its coor- 
dinates do not satisfy any algebraic equation with coefficients in Q besides those that are 
satisfied by every point on V. 

A point x on a semi-algebraic set S is locally generic if for small enough e, x is generic in 
S(lB t (x). 

Remark 2.4. A semi- algebraic set S will have no generic points if its Zariski closure over Q 
is reducible over Q. In this case all points in each irreducible component will satisfy some 
specific equations not satisfied everywhere over 5*. But even in this case, almost every point 
in the set will be generic within its own component. Such points will still be locally generic 
in S. 

We will need a few elementary lemmas on generic points in semi-algebraic sets. 

Lemma 2.5. // X and Y are both semi-algebraic sets defined over Q, with X C Y and 
X dense in Y (in the Euclidean topology), then X contains all of the locally generic points 
ofY. 

Proof. Due to the density assumption, Y\X must be a semi-algebraic set defined over Q 
with dimension less than that of Y. The Zariski closure of a semi-algebraic set maintains 
its dimension, thus all points in Y\X must satisfy some algebraic equation that is non-zero 
over Y and thus these points must be non-generic. To see that a point y G Y \ X is not 
locally generic either, apply this argument in an e-neighborhood of y. (That is, apply the 
argument to Y fl B e (y) and X fl B e (y).) □ 
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Lemma 2.6. Let V be a semi- algebraic set, f be an algebraic map from V to a range space X, 
and W be an semi-algebraic set contained in the image of f , with V , W , and f all defined 
over Q. If xq G V is generic and f(xo) G W, then f(xo) is generic in W. 

Proof. Let <j) be any algebraic function on X with rational coefficients so that <j)(f(xo)) = 0. 
Then <fi o / vanishes at Xo, which is generic in V, so <f> o f vanishes identically on V. This 
implies that <fi vanishes on W. Thus any algebraic function defined over Q that vanishes at 
f(xo) must vanish on all of W. This proves that f(xo) is generic in W. □ 

3. The geometry of PSD stresses 
We now turn to the main construction in our proof. 
3.1. The measurement set. 

Definition 3.1. The d-dimensional measurement set M^(T) of a graph T is defined to be the 
image of C d (V) under the map £r- These are nested by MdiV) C Md+i(T) and eventually 
stabilize at M„_i(T), also called the absolute measurement set M(T). 

Since Md(T) is the image of C d (V) under an algebraic map, by Lemma 2.6, if p is generic 
in C d {V) then £r{p) is generic in Md(T) (which must be irreducible). 

Lemma 3.2. The set M(T) is convex. 

Proof. The squared edge lengths of a framework p are computed by summing the squared 
edge lengths of each coordinate projection of p, so MdiT) is the <i-fold Minkowski sum of 
Mi(T) with itself. Since M\{T) is invariant under scaling by positive reals, Md(T) can also 
be described as an iterated chord variety, and in particular M(T) is the convex hull. □ 

A particular case of interest is M(A V ), the absolute measurement set of the complete 
graph on the vertices. It is a cone in W^ 2 \ 

Lemma 3.3. The measurement set Md(A v ) is isomorphic (as a subset of the linear space 
Rw ) to the set of PSD (v — 1) x (y — 1) matrices of rank at most d. 

Proof. There is a standard linear map from M(A V ) to the convex cone of v — 1 by v — 1 
positive semidefinite matrices, which we recall for the reader. (See, e.g., [Sch35,Gow85].) For 
a point x G Md(A v ), by the universal rigidity of the simplex, up to <i-dimensional Euclidean 
isometry there is a unique p G C d (V) with x = £a(p)- If we additionally constrain the last 
vertex to be at the origin, then p is unique up to an orthonormal transform. Think of such 
a constrained p as a (t> — 1) x d matrix, denoted g. The map sends x to the matrix gg t , 
which has rank at most d. (The rank of gg l is less than d if p has an affine span of dimension 
less than d). It is easy to check that gg l does not change if we change g by an orthonormal 
transform, and that in fact the entries of gg t are linear combinations of the squared lengths 
of edges, so the map is well-defined and linear. □ 

Definition 3.4. Let K be a non-empty convex set in an affine space. The dimension of K is 
the dimension of the smallest affine subspace containing K. A point x G K is in the relative 
interior Int(-fT) of K if there is a neighborhood U of x in E e so that U H K = M. dimK . (It is 
easy to see that the relative interior of K is always non-empty.) 

The face F(x) of a point x G K is the (unique) largest convex set contained in K so that 
x G Int(F(x)). The faces of K are the sets F(x) as x ranges over K. (In Section 4.1 we will 
see more characterizations of F(x).) Let f(x) be the dimension of F(x). 
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Lemma 3.5. Let x = £ A (p) with p G C d (V). Then F(x) = £ A (AS d (p)) = £ A {GL d (p)). 
When p is proper, f(x) = r 2 ^ 1 

Here GL d is the group of d- dimensional linear transforms, and Aff^ is the corresponding 
affine group. 

Proof. By Lemma 3.3, this reduces to understanding the faces of the cone of v — 1 by v — 1 
PSD matrices, which is well understood (see, e.g., [PatOO, Theorem A.2]). In particular, for 
q a (v — 1) X d matrix, the points in F(gg l ) in the semidefmite cone are those matrices of 
the form gL t Lg t for some dx d matrix L. But each such gL t Lg t maps under our isometry to 
£ A (a) where a G C d (V) is obtained from p by a d- dimensional affine transform, with linear 
part described by L. When p has an affine span of dimension d, g is of rank d. This implies 
that F(gg t ) has dimension l^ 1 )' as we can see by computing the dimension of GL d /O d , the 
linear transforms modulo the orthonormal transforms. □ 

3.2. Dual vectors and PSD stresses. 

Definition 3.6. Let us index each of the coordinates of with an integer pair ij with 
1 < 2 < j < f . Given 0, a functional in the dual space (mW^ , define the v-by-v matrix 
M(0) as follows: for i ^ j, M(0)^- = M(0)^ := 0y and M(0) il := - £^ M(0)y 

Lemma 3.7. Let be a dual vector in (W^^j and let f2 := M(0) be the corresponding 
matrix. LetpEC d (V). Then 

(3-8) (0,Mp)> = t^X>W- 

k=l 

Here we use the notation p k for vector in M. v describing the component of p in the /c'th 
coordinate direction of E d . 

Proof. Both sides are equal to J2k^2i,j\i<ji.Pi ~ PjY'Pij- D 

Definition 3.9. For K C W 1 a convex set and a functional in the dual space (M n )*, we 
say that is tangent to K at x 6 K if for all y G K, ((f), y) > while (0, x) = 0. 

Lemma 3.10. Let (j) be a dual vector in (r( 2 )^ that is tangent to M(A V ) at £a(p) for some 
p G C(V). Let Q := M(0) be the corresponding matrix. Then Q is a PSD equilibrium stress 
matrix for (p, A„) . 

Proof. Conditions (1) and (3) in Definition 1.8 are automatic by definition of Vt, and condi- 
tion (2) is vacuous for the complete graph A„. It remains to check condition (4) and that Q 
is PSD. 

By Lemma 3.7, (0, •) can be evaluated as a quadratic form using the matrix Q. Since (0, •) 
is not negative anywhere on M(A V ), Q must be PSD. Because Q is positive semidefinite and 
J2k=i{p k y^P k = 0, it must also be true that for each k, (p k ) l VLp k = and so Qp k = 0, which 
is the last necessary condition to show Q is an equilibrium stress matrix. □ 

Lemma 3.11. In the setting of Lemma 3.10, suppose furthermore that <fr is only tangent to 
points in F(£ A (p)), and suppose the affine span of p has dimension d. 

Then fl is a PSD equilibrium stress matrix for (p, A v ) with rank v — d — 1. 
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Proof. From Lemma 3.10, Q is a PSD equilibrium stress matrix for (p, A v ). By assumption, 
its kernel is spanned by the coordinates of frameworks that map under £a to F From 
Lemma 3.5, such frameworks are <i-dimensional affine transforms of p. Thus the kernel is 
spanned by the d coordinates of p and the all-ones vector. Since the kernel has dimension 
d + 1 , Q has rank v — d — 1 . □ 

Lemma 3.12. In the setting of Lemma 3.11, suppose furthermore that <pij = for {i,j} ^ 
£(r) for a graph T. 

Then Q is a positive semidefinite equilibrium stress matrix for (p, T) with rank v — d — 1. 

Proof. Q must have zeros in all coordinates corresponding to non-edges in T, thus it will be 
an equilibrium stress matrix for T as well as for A„. The rest follows from Lemma 3.11. □ 

To be able to use Lemma 3.12, we want to find a dual vector that is only tangent to 
points in F(£^(p)) C M(A V ). Additionally we need that <pij = for {i, j} G" S. As we will 
see, this will correspond to finding a dual vector in (IR e )* that is only tangent to points in 
F(£r(p)) C M(r). The proper language for this kind of condition is the notion of convex 
extreme and exposed points, which we turn to next. 

4. Convex sets and projection 

We will prove our main theorems in the more general setting of closed convex cones. Our 
goal is the following general theorem about projections of convex sets, which will quickly 
imply Theorem 1. (See Definitions 4.1, 4.6, and 4.11 for the terms used.) 

Theorem 2. Let K be a closed line-free convex semi-algebraic set in M™ 1 , and ir : W 71 — > M n 
a projection, both defined over Q. Suppose x is locally generic in extk(K) and universally 
rigid under it. Then tt(x) is k-exposed. 

4.1. Extreme points. 

Definition 4.1. Let K be a non-empty, convex set. A point x G K is k-extreme if f(x) < k. 
(Recall from Definition 3.4 that f{x) is the dimension of the face of K containing x.) Let 
extk(K) be the set of fc-extreme points in K. It is easy to see that extk(K) C extk+i(K) 
and if K is closed, then extfe(i^) is closed. 

We will also use the following elementary propositions (see for example the exercises in 
chapter 2.4 of [Gru03].) 

Proposition 4.2. For K a convex set and x G K , the following statements are equivalent: 

• f(x) < k, and 

• x is not in the relative interior of any non- degenerate (k + 1) -simplex contained in K. 

Proposition 4.3. For K a convex set and x G K, the face F(x) is the set of points z G K 
so that there is an z' G K with x in the relative interior of the segment [z f , z] . 

Remark 4.4. One special case of Proposition 4.3 is when z = x, in which case we also take 
z' = x and the segment consists of a single point. 

Corollary 4.5. For K a convex set, the faces of K are the convex subsets F of K such that 
every line segment in K with a relative interior point in F has both endpoints in F . 
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Figure 2. Examples of extreme and exposed points. For the racetrack shape 
illustrated (bounded by semi-circles and line segments), a is 0-extreme and 
0-exposed, b is 1-extreme and 1-exposed, c is 2-extreme and 2-exposed, and d 
is 0-extreme and 1-exposed. 

4.2. Exposed points. 

Definition 4.6. A point x £ K is k- exposed if there is a closed half-space H containing x 
so that dimH n K < k. Let exp k (K) be the set of fc-exposed points in K. 

See Figure 2 for some examples of fc-extreme and fc-exposed points, including a case where 
they differ. If x is fc-exposed, it is also /c-extreme, as any simplex containing x in its relative 
interior is contained in any supporting hyperplane. 

The following theorem is a crucial tool in our proof. 

Theorem 4.7 (Asplund [Asp63]). For K a closed convex set, the k- exposed points are dense 
in the k-extreme points. 

Remark 4.8. Asplund only states and proves Theorem 4.7 for compact convex sets. The 
result follows for any closed convex set K as follows. For x £ ext k (K), fix R > and 
consider the compact convex set K' :— K H Br{x). Then Asplund's theorem says we can 
find /c-exposed points z in K' that are arbitrarily close to x. If z is sufficiently close to x 
(say, within R/2 of x), it will also be /c-exposed in K, as desired. 

Theorem 4.7 is an extension of a theorem of Straszewicz [Str35], who proved it in the 
case k = 0. There are several improvements of Theorem 4.7, replacing "dense" with various 
stronger assertions. In our applications we will consider semi-algebraic sets K defined over Q. 
Note that for any such set K, extk(K) and exp k (K) are also semi-algebraic over Q, as the 
/c-extreme and /c-exposed conditions can be phrased algebraically. 

Corollary 4.9. Let K be a closed, convex, semi- algebraic set. If x is locally generic in 
extk(K), then x is k-exposed. 

Proof. This follows from an application of Lemma 2.5 to the inclusion exp k (K) C ext k {K) 
in a small neighborhood around x. □ 

4.3. Projection. The two convex sets M(A V ) and M(T) are related by a projection from 

r( 2 ) to 1R 6 . We will consider this situation more generally. 

Throughout this section, fix two Euclidean spaces, E m and E n , with a surjective projection 
map 7r : E m — > E n and a closed convex set K C E m . We will work with the image w(K) C E n . 
Let 7Tk be the restriction of ir to K. 

Remark 4.10. In general, tt(K) need not be closed, even if K is. We continue to work with 
n(K) and its "faces" as defined in Definition 3.4, even if it is not closed. 
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Definition 4.11. We say that x £ K is universally rigid (UR) under tt if ir K (tt(x)) is the 
single point x. 

Note that for y £ 7r(fT), 7r^ (y) is always a convex set in its own right. 

Lemma 4.12. Let K be a convex set and it a projection. For any x £ K, ir(F(x)) C 
F(<k(x)). 

Proof. The set tt(F(x)) is a convex set containing x in its relative interior, so by definition 
is contained in F(tt(x)). □ 

A point that is 0-extreme in a convex set is called a vertex of the set. A convex set is line 
free if it contains no complete affine line. Recall that a non-empty, closed, line-free convex 
set has a vertex (see e.g., [Grii03, 2.4.6]). Our convex sets of interest, M(Ka) and M(T), 
are closed and also automatically line free, as the length- squared takes only positive values 
and any line contains both positive and negative values in at least one coordinate. 

Proposition 4.13. Let K be a convex set and n a projection. For any y £ tt(K) and x a 
vertex of 7r]^(y), F(x) maps injectively under tik into F(y). In particular, f(x) < f(y). If 
K is closed and line free, then for every y £ tt(K) there is a x £ ir~ 1 (y) so that f(x) < f(y)- 

Proof. Let A(x) be the smallest affine subspace containing F(x). Suppose A(x) contained 
a direction vector v in the kernel of n. Then, since x is in the relative interior of F(x) for 
small enough e, we would have the segment [x + ev, x — ev] contained in F(x) and also in 
7i] ( 1 (y). This would contradict the assumption that x is a vertex of ^^(y)- Thus F(x) maps 
injectively under tik- 

From Lemma 4.12 we see that f(x) < f{y). For the last part of the lemma, observe that 
if K is closed and line free, so is (y), and so ^^(y) has a vertex, which is the desired 
point x by the first part. □ 

Proposition 4.14. Let K be a convex set and ir a projection. For any y £ n(K) and x in 
the relative interior ofn I< 1 (y), F(x) = ixJ^(F(y)). In particular, f(x) > f(y). 

Proof. By Lemma 4.12, we already know that F(x) C tt^ 1 (F (y)) , so we just need to show 
the other inclusion. Pick any y 1 £ F{y) and X\ £ We must show that x\ £ F{x). 

Since y\ £ F(y), by Proposition 4.3 there is an y 2 £ n(K) so that y £ IntQyi, 2/2])- Let 
£2 £ t^k (2/2)- Let x' £ K be the (unique) point of intersection of [xi,X2] with ^^(y). Since 
x £ Int^^ 1 (y)), there is a point x" £ ^^(y) with x £ Int([x', x"]). But then x is in the 
interior of the simplex [xi,X2,x"], which implies that x\ £ F(x). □ 

In order to be able to apply Asplund's theorem at n(x) we need n(x) to be locally generic 
in extfc(7r(J^)). The following lemma will assist us. 

Lemma 4.15. Let K be a closed, convex set in IR m and let %k '■ K — >■ lR n be a projection. Let 
x £ K be universally rigid under ir. Then for any e there is a 6 > so that tx]^{Bs{tx{x))) C 
B e (x) n K. 

Intuitively, the only points in K that map close to ir(x) are close to x. 

Proof. Suppose not. Then there is an e and a sequence of points Xi with the following 
property: 7r(xj) approach ir(x), while d(x,X{) > e. Let x\ be the point on the interval [x,Xi] 
that is a distance exactly e from x. By convexity x\ £ K. Then the x\ are in a bounded and 
closed set and therefore have an accumulation point x' £ K, which will also have distance e 
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from x. By the linearity of tt the 7r(a^) also approach n(x). So by continuity ir(x') = ir(x). 
This contradicts the universal rigidity of x. □ 

And now we can prove the following. 



Lemma 4.16. Let K be a closed line-free convex semi-algebraic set in M. m and tt : M. m — > W n 
a projection, both defined over Q. Suppose x is locally generic in extfc(A) and universally 
rigid under it. Then ir(x) is locally generic in extfc(7r(A)). 

Proof. By Proposition 4.13, extfc(7r(A)) C 7r(extfc(A)). By Proposition 4.14 and the universal 
rigidity of x, we have that tt(x) is in extfc(7r(A')). 

Let V := extfc(A) fl B e (x) and W := ext fc (7r(A)) fl Bs{n(x)). For sufficiently small e, by 
local genericity, x is generic in V. Meanwhile, W C 7r(extfc(A)), and, from Lemma 4.15 , for 
small enough 5, we have ^^(W) C B e (x) and thus W C 7t(V). 

Thus from Lemma 2.6, we have tt(x) generic in W. Thus ir(x) is locally generic in 
ext fe (7r(A)). □ 

Remark 4.17. In our special case where A is the cone of positive semidefinite matrices, x is 
in fact generic in extfc(A) (which is irreducible), and thus the second paragraph in the above 
proof (and in turn Lemma 4.15) are not needed. 

Finally, in order to be able to apply Asplund's theorem to tt(K) we need tt(K) to be a 
closed set. In the graph embedding case, we can use the fact that it is a proper map whenever 
T is connected, and thus ir(K) must be closed. In the setting of a general A and n we can 
argue closedness using standard techniques. 

Definition 4.18. A direction of recession v for a convex set A is a vector such that for every 
(equivalently, any) x G A and every A > we have x + Xv e A (see e.g. [Roc70, Chapter 
8]). 

A sufficient condition for closedness is given by the following theorem [Roc70, Theorem 
9.1]. 

Theorem 4.19. Let A be a line free closed convex set and ir a projection. If A does not 
have a direction of recession v with ir(v) = 0, then then 7r(A) is closed. 

Corollary 4.20. Let A be a line free closed convex set and it a projection. If there is a 
point x G A that is universally rigid under tt, then ir(K) is closed. 

Proof. Suppose there was a direction of recession v for A with tt(v) = 0. Then tt(x+v ) = ti(v ) 
with x + v in A. This would contradict the universal rigidity of x. □ 

Putting this all together, we can deduce Theorem 2. 

Proof of Theorem 2. From Corollary 4.20, 7r(A) is closed. From Lemma 4.16, ir(x) is locally 
generic in extfc(-7r(A)). Thus by Corollary 4.9, ir(x) is ^-exposed. □ 

4.4. Proof of main theorem. We are now in position to complete the proof of Theorem 1. 
Recall that we have a universally rigid framework (p, T) with p generic in C d (V). Let A be 
M(A V ), and vr be its projection to M(T), k be and x G M d (A v ) C A be £ A (p). The 

assumption that p is generic in C d (V) implies that x is generic in Md(A v ) = extfc(A). 

From Theorem 2, ir(x) is /c-exposed. So there must be a halfspace H in R e whose inter- 
section with 7r(A) is exactly F(tt(x)). 
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The preimage n~ 1 (H) is a halfspace in Rv 2 ) whose intersection with K is n] < 1 (F(Tr(x))), 
which by Proposition 4.14 is F{x). 

When v > d+1, 7r~ 1 (if) must be tangent to K. Since K is a cone, the boundary of ir~ l (H) 

must pass through the origin. Thus tt~ 1 (H) is represented by a dual vector G (iR^ 2 )) , in 

the sense that 7T _1 (if) = { x G rGO <f>(x) > }. 

Since we started with a halfspace in M e , 0y = for ^ £• From Lemma 3.12, then, 

f2 := M(0) must be a positive semi-definite equilibrium stress matrix of (p, T) with rank 
v — d — 1. 

5. Relation to semidefinite programming 

Theorem 1 can be interpreted as a strict complementarity statement about a particular 
semidefinite program (SDP) and its associated dual program, as we will now explain. For 
a general survey on semidefinite programming see Vandenberghe and Boyd [VB96]. Our 
notation is similar to that of Pataki and Tungel [PT01]. For a related discussion on universal 
rigidity, semidefinite programming and complementarity, see, for instance, the paper by So 
and Ye [SY07]. 

5.1. Semidefinite programming. We first recall the basic definitions of semidefinite pro- 
gramming. Let S n be the linear space of all symmetric n-bj-n matrices and S 1 ™ C S n 
be the cone of positive semidefinite matrices. A semidefinite program is given by a triple 
(L,b,/3) where L C S n is a linear subspace, b G S n is a matrix, and (3 G (S n )* is a 
dual to a matrix (which we can also view as a matrix, using the natural inner product 
(A, B) = tr(A t B) = AijBij). This describes a primal constrained optimization problem: 

inf{(/3,x) \ x e(L + b)nSl} 

That is, we optimize a linear functional over all symmetric positive semidefinite matrices x, 
subject to the constraint that x lies in a chosen affine space. 

Associated with every primal SDP is the (Lagrangian) dual program 

inf{(fi,6) | Q G {L L + (5) n (SI)*} 

Since the cone of PSD matrices is self-dual, elements of (S 1 ")* also correspond to PSD 
matrices, so this dual program can be thought of as a semidefinite program itself. 

Definition 5.1. We say a pair of primal/dual points (x, Q) with and Q G (S 1 ™)* are 

a complementary pair for the program (L, b, /3) if both are feasible points for the respective 
programs and (Q,x) = 0. 

A calculation shows that for any feasible primal/dual pair of points (x, Q), we have ((3, b) — 
(Q,b) < ((3,x) and that complementarity implies {/3,b) — {fl,b) = ((3,x). 
Thus, writing the dual problem in the alternative form 

su P {(/3, b) - <n, b) | n g (l x + 0) n (£?)*} 

n 

we see that a complementary pair must represent optimal solutions to the primal and dual 
SDPs respectively. 

Definition 5.2. We say that an SDP problem (L, b, (3) is gap free if it has complementary 
pairs of solutions. 
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Gap freedom is typically seen as a fairly mild constraint on an SDP problem. 
For every complementary pair (x, fi) we have [AH097] 

(5.3) rank(x) + rank(f2) < n. 

Definition 5.4. A complementary pair (x, Q) is said to be strictly complementary if we have 
rank(x) + rank(f2) = n. 

There is a rich theory on when SDPs have strictly complementary solution pairs. (See 
[AHO97,PT01,NRS10]; see also Section 5.9.) In particular, certain SDP algorithms converge 
more quickly on problems that satisfy strict complementarity [LSZ98]. The linear program- 
ming counterpart of strict complementarity is strict complementary slackness [DT03]. In 
the linear case, strict complementarity is always achieved. By contrast, for SDP problems it 
depends on the particular problem. 

5.2. Graph embedding as an SDP. In the context of graph embedding, suppose we are 
looking for a configuration of unconstrained dimension, p £ C V (V)/ Eucl(u), such that the 
framework of V is constrained to have a set of squared edge lengths d?j for ij £ £(r). It is well 
known that we can set up this graph embedding problem as a semidefinite program [LLR95], 
as we will now review. 

To a configuration p we can associate a v x v Gram matrix x, with Xij = p(i) ■ p{j). The 
Gram matrix of p is unchanged by elements of 0(v) (although it does change when p is 
translated). The matrix rank, r, of such an x is the dimension of the linear span of the 
associated configuration p. This will be typically be d + 1; one greater than the affine span, 
d, of the framework. 

In the graph embedding SDP, we set n := v and take the Gram matrix x as our unknown. 
The distance constraint at an edge ij can be expressed as the linear constraint xa+Xjj—2xij = 
{dij) 2 . The collection of matrices satisfying these constraints for all edges forms an affine 
space; we choose L and b such that L + b is this affine space. In our context, we are only 
interested in feasibility, and thus have no objective function, so we set f3 := 0. Note that 
semidefinite programs do not allow us to explicitly constrain the rank of the solution (which 
is related to the dimension of the configuration). 

Let us now look at the dual program to our graph embedding SDP problem. In the primal 
problem, the linear space L corresponds to symmetric matrices x with xu + Xjj — 2xij = 
for all {i, j} £ £(r). Thus the space L- 1 , when represented as matrices, is spanned by a basis 
elements = —en + e^ + — Cjj for {i,j} £ £(r), where e^- is the elementary matrix 
with a 1 in the ij entry and elsewhere. The space L x is therefore those matrices Q with 
row-sums of zero and with Q.^ = for all % ^ j, {i,j} ^ £(r). That is, L x is the space of all 
(not necessarily equilibrium) stress matrices. We also have /3 = 0, so we optimize over these 
matrices. 

Minimizing (Q, b) imposes the further constraint that the solution be an equilibrium stress 
matrix. This can be seen as follows. Since ,5 = 0, the all-zeros matrix is dual feasible. 
Whenever there is some primal feasible x corresponding to a configuration p, the SDP must 
be gap free as (Q,x) = 0, and the optimal Q are the feasible Q with (Q,x) = 0. For any 
\& £ (S 1 ™)* we have {$,x) = \ Y^k=iip k ) t ^P k ■ (Here we use the notation p k for vector in 
W 1 describing the component of p in the fc'th coordinate direction of lR r .) And thus, as in 
Lemma 3.10, the feasible Q with (Q,x) = are the PSD equilibrium stress matrices for 

(p,r). 

In this language we can restate Theorem 1 as follows. 
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Proposition 5.5. Let p be generic in C d (V), let (p,T) be universally rigid in E , and let 
(L,b,0) be the associated SDP using the graph V and the distances of the edges in (p, T). 
Then (L, b, 0) has a strictly complementary solution pair (x,Q). 

Proof. Simply let x be the Gram matrix corresponding to a translation of p such that its affine 
span does not include the origin. Thus x will have rank d + 1. Let Q be a PSD equilibrium 
stress matrix of rank v — d — 1, which exists by Theorem 1. Then rankx + rankfi = v. □ 

5.3. SDP feasibility. Our discussion of universal rigidity and equilibrium stress matrices 
carries over directly to any SDP feasibility problem (L, b, 0) (i.e., where (3 = 0). To see 
this we first set up some notation. Pick an integer r > and let k := ( r ^ X )- Then, as in 
Lemma 3.5, extfc(S'") is the set of the PSD matrices of rank less than or equal to r, which 
we will denote S™ r . Let n be projection from S n to S n /L, and ir + its restriction to S*™. A 
point x G S 1 } is a solution to the feasibility problem iff tt(x) = ir(b), or x G 7r^ 1 (7r(6)). 

A point x G S™ is universally rigid under n iff 7r^ 1 (7r(x)) is a single point: A universally 
rigid feasible point is the unique feasible solution to (L, 6, 0). As we vary b with fixed L, we 
will see as solutions all points in £?; it may happen that for some b there is a unique, generic 
feasible solution. 

Proposition 5.6. Let (L, b, 0) be an SDP feasibility problem, where L has rational coeffi- 
cients. Suppose there is a unique feasible solution x which is generic in S™ r . Then there 
exists an optimal dual solution Q such that (x, Q) is a strictly complementary pair. 

Proof. Since x is generic in 5*",, = extfc(S'"), by Theorem 2 ir(x) is /c-exposed. Thus there is 
a PSD matrix Q in (S n /L)* = L 1 - that is tangent to ir(S+) at 7r(x), with contact F(tt(x)). 
Since Q G it is a feasible point of the dual SDP. Since Q is tangent to tt(S^) at 7r(x), it 
is also tangent to S™ at x; thus (Q, x) = and O is a complementary and therefore optimal 
dual solution. 

The tangency of f2 to 5*" at x has contact 7r^ 1 (F(7r(x)) which, by Proposition 4.14, is F(x). 
Let p be an r-dimensional configuration of n points with Gram matrix x. From the facial 
stucture of the PSD cone we see that F(x) consists only of Gram matrices of r-dimensional 
linear transforms of p. We also have the relation (O, x) = ~ J2k=i{p k ) t &'P k = 0- Thus we see 
that Q must have a kernel of dimension r and have rank n — r. □ 

5.4. SDP optimization. We can apply the approach above to a more general SDP problem 
(L,b,(3), where /3 ^ 0. Specifically, we prove the following. 

Theorem 3. Let (L, b, (3) be a gap-free SDP problem, with L and (3 rational. Suppose 
there is a unique optimal solution x which is generic in S^L.. Then there exists a strictly 
complementary pair (x, Q) . 

We will prove the theorem by reducing to the (3 = case. For any SDP optimization 
problem (L, b, (3) with an optimal solution x, there is an associated SDP feasibility problem 
(17, b\ 0), where LI + b' is {y G L + b | y) = x)}. In particular, L' is L R ker (3 (which 
is still rational) and the dual space L /J - is L L + ((3). Then feasible solutions to (L',b',0) 
correspond to optimal solutions to (L,b,(3). 

Lemma 5.7. Let (L,b,(3) be an SDP problem, with L and (3 rational. Suppose there is a 
unique optimal solution x which is generic in 5™ r . Then there exists a PSD Q such that 
rankf2 = n — r, with Q G L L + ((3) and (Q,x) = 0. 
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Proof. Apply Proposition 5.6 to the feasibility problem (L', b', 0) constructed above. This 
gives a Q G L = L 1 - + ((3) that is strictly complementary to x. □ 

Note that we have not yet proved Theorem 3, as Lemma 5.7 gives fl in the linear space 
L + (j3) rather than the desired affine space L + j3. 

Proof of Theorem 3. Let Qi be the PSD matrix given by Lemma 5.7. From the assumption 
of gap freedom there exists a PSD Q 2 with Q 2 £ L 1 - + (3 and (Q 2 ,x) = 0. Thus for any 
positive scalars Ai and X 2 , the matrix Q := Aif^ + \ 2 Q 2 is PSD, has rank no less than n — r, 
and has (Q, x) = 0. By adjusting Ai and \ 2 we can achieve Q G L 1 - + j3. By Equation (5.3), 
we in fact have rank Q = n — r. □ 

We note though that for a given (L, (3) and rank r, as we vary b, there may be no unique 
solutions with rank r. Even if there are such x, they may all be non-generic. In such cases, 
Theorem 3 is vacuous, as there are no choices of (L, b, (3) satisfying the hypotheses. 

5.5. Genericity. In the context of SDP optimization, the hypothesis in Theorem 3 that L 
and f3 be rational may seem a little unnatural. This hypothesis can be relaxed if we work 
with points that are generic over a field that is larger than Q. 

Definition 5.8. Let k be a field containing Q and contained in IR. A semi- algebraic set S C 
W 1 is defined over k if there is a set of equalities and inequalities defining S with coefficients 
in k. If S is defined over k, a point x G S is generic over k if the coordinates of x do 
not satisfy any algebraic equation with coefficients in k beyond those that are satisfied by 
every point in S. Similarly, x is locally generic over k if for small enough e, x is generic in 
S(lB e (x). 

A defining field of a semi-algebraic set S, written Q[S], is any field over which it is defined. 
(There is a unique smallest field for algebraic sets by a result of Weil [Wei46, Corollary IV. 3]. 
For our purposes, we can use any field over which S is defined.) Similarly, if / : X — > Y is a 
map between semi-algebraic sets, then Q[f] is a field over which it is defined (or, equivalently, 
a field over which the graph of / is defined). 

For instance, if S is a single point x, we can take Q[S] to be the same as Q[x], the smallest 
field containing all the coordinates of x. Also, if x is generic over Q, then y is generic over 
Q[x] iff the pair (x, y) is generic over Q. 

With this definition, Theorems 2 and 3 can be improved to allow non-rational sets and 
projections. 

Theorem 2'. Let K be a closed line-free convex semi-algebraic set in MJ 71 , and ir : IR m — > R n 
a projection. Suppose x is locally generic over Q[K,ir] in extk(K) and universally rigid 
under ir. Then ir(x) is k-exposed. 

Theorem 3'. Let (L, b, (3) be a gap-free SDP problem. Suppose there is a unique optimal 
solution x which is generic over Q[L,(3] in S™ r . Then there exists a strictly complementary 
pair (x, fi) . 

The proofs follow exactly the proofs of the versions given earlier. For instance, in Theo- 
rem 3', we apply Theorem 2' to the cone 5+ (which is defined over Q) and the projection 
onto S n /(L + (/?)), which is defined over whatever field is needed to define L and (3. 
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5.6. Relation to previous results. A fundamental result of [AHO97,PT01] on strict com- 
plementarity can be summarized as follows. 

Theorem 5.9. Suppose (L,b,f3) is a generic SDP program that is gap free. Then (L,b,/3) 
admits a strictly complementary pair of solutions. 

This result is neither stronger nor weaker than our Theorem 3. In particular Theorem 5.9 
requires that all parameters be generic. In contrast, Theorem 3 does not assume genericity of 
any of the parameters but rather assumes genericity of the solution within its rank. Indeed, in 
our application to rigidity of graphs, ,5 = 0, which is obviously not generic. (The parameter 
b is also not usually generic.) In fact, there are very few situations where both theorems can 
apply, as the following proposition shows. 

Proposition 5.10. In an SDP problem (L,b,(3) with dimL = D and primal solution x of 
rank r, if Theorem 5.9 applies, then 

1 qn — —I— 1 \ 

< D. 



> D. 




On the other hand, if Theorem 3 applies, then 

(n — r + 1 
I 2 



Proof. For a given (L,b,(3), if b is generic over Q[L], then r satisfies the first inequality, as 
described in [AH097, Theorem 12] and [NRS10, Proposition 5]. In particular, for generic b, 
the intersection of L + b with S?t r must be transversal, which from a dimension count gives 
the first inequality. This takes care of the first part of the proposition. (If /3 is generic over 
Q[L], as in Theorem 5.9, there is also an upper bound on r: ( r ^ 1 ) — ("2^) — ^' we ^° 
not need this.) 

For the second part, recall that if Theorem 3 applies, there is a point x, generic over 
Q[L, (3] in S™ r which is the unique solution to the SDP problem (L, b, (3). Recall that there is 
an associated feasibility problem (L', b', 0). Let n : S n — > S n /L' be the associated projection 
and 7r + its restriction to S 1 ™. Since x is unique (in both (L,b,(3) and (L',6',0)), x is the 
only point in S 1 " mapping to ir(x), and tt(x) G d7r(S™). Moreover, since x is generic in S™ r , 
there is an open neighborhood U of a; in S™ r with these properties. In particular, U maps 
injectively by n to dn(S$) C S n / L' . 

3%, has dimension - ( n ~2 +1 ). On the other hand, dim(5 n /L') < ( n + x ) -D + l 

(with equality iff (3 & or equivalently V ^ L), so dim(97r(S'™)) < ( n + 1 ) - D. If there is 



a smooth injection from S™ r to 9^(5™), we must have ( n ^ +1 J > D, as desired. □ 

5.7. Complexity of universal rigidity. Assuming (p,T) is not at a conic at infinity and 
is translated so that its afline span does not include 0, (p, T) is UR iff there is no higher 
rank solution to the SDP than p [Alf07a]. Thus, to test for universal rigidity, the main step 
is to test if there is a feasible solution with rank higher than that of a known input feasible 
solution p (given say as integers). Numerically speaking, SDP optimization algorithms that 
use interior point methods produce approximate solutions of the highest possible rank [GY93, 
DKRT97] and so in practice one could try such a method to produce a guess about the 
universal rigidity of (p, T). 
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The complexity of getting a definitive answer is a trickier question, even assuming one 
could reduce the UR question to one of "yes-no" SDP feasibility. In particular, an approx- 
imate solution to an SDP optimization problem can be found in polynomial time but the 
complexity of the SDP feasibility problem remains unknown, even with strict complemen- 
tarity. See [Ram97] for formal details. 

Theorem 1 tells us that for generic inputs, the UR question can be answered by finding 
the highest rank dual optimal solution. (A framework with integer coordinates will not be 
generic; however, if the integers are large enough we are likely to avoid all special behavior, as 
in [GHT10, Section 5].) This does not appear to be any help in determining the complexity 
of testing algorithmically for universal rigidity. 
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